library=
{
	{ "H+",   1 },
	{ "HO-", -1 },
	{ "HA",   0 },
	{ "A-",  -1 },
	{ "Na+",  1 }
};

chemsys =
{
	{ "water",   10^(-14),  { 1, "H+" }, { 1, "HO-" } },
	{ "vinegar", 10^(-4.8), {-1, "HA" }, { 1, "H+" }, { 1, "A-" } }
};


C0 = 0.01;
ini =
{
	{ C0,  { 1, "HA" }, { 1, "A-"} }, -- initial concentration
	{ 0.00, { 1, "Na+" } }, -- no counter ion
	{ 0.00, { 1, "H+" }, {-1,"HO-"}, { -1, "A-" }, { 1, "Na+" } } -- electroneutrality to start with 
}

Volume = 20.00; -- mL
Cb     =  0.01; -- M
N0     = C0 * Volume;
Veq    = N0/Cb;
Vb     = 2*Veq;
